Crystallographic Studies of KRAS in Complex with Small Molecules and RAS-Binding Proteins.

RAS proteins play a vital role in regulating downstream signaling and essential cellular processes, positioning them as key players in normal cellular physiology and disease development. Among the various isoforms of RAS, KRAS stands out as one of the most frequently mutated genes in human cancer. The prevalence of RAS mutations in cancer often involves single amino acid substitutions at codons 12, 13, or 61. These mutations disrupt the RAS protein's inherent ability to transition between its active and inactive states, resulting in a constant activation signal and driving uncontrolled cell growth. Crystallization and structural analysis of KRAS with inhibitors and RAS-binding proteins play a pivotal role in unraveling the structural and mechanistic details of KRAS function, aiding in drug discovery efforts, and advancing our understanding of KRAS-driven diseases. Here, we present our experimental methodology for crystallizing KRAS in the presence of covalent or non-covalent small molecules and proteins acting as effectors or regulators of RAS. We detail the techniques for successful crystallization and the subsequent optimization of crystallization conditions. The resulting crystals and their structures will provide valuable insights into the key interactions between KRAS and its partner proteins or potential inhibitors, offering a foundation for developing targeted therapies that are more potent and selective against KRAS-driven cancers.

Structural dynamics of RAF1-HSP90-CDC37 and HSP90 complexes reveal asymmetric client interactions and key structural elements.

RAF kinases are integral to the RAS-MAPK signaling pathway, and proper RAF1 folding relies on its interaction with the chaperone HSP90 and the cochaperone CDC37. Understanding the intricate molecular interactions governing RAF1 folding is crucial for comprehending this process. Here, we present a cryo-EM structure of the closed-state RAF1-HSP90-CDC37 complex, where the C-lobe of the RAF1 kinase domain binds to one side of the HSP90 dimer, and an unfolded N-lobe segment of the RAF1 kinase domain threads through the center of the HSP90 dimer. CDC37 binds to the kinase C-lobe, mimicking the N-lobe with its HxNI motif. We also describe structures of HSP90 dimers without RAF1 and CDC37, displaying only N-terminal and middle domains, which we term the semi-open state. Employing 1 µs atomistic simulations, energetic decomposition, and comparative structural analysis, we elucidate the dynamics and interactions within these complexes. Our quantitative analysis reveals that CDC37 bridges the HSP90-RAF1 interaction, RAF1 binds HSP90 asymmetrically, and that HSP90 structural elements engage RAF1's unfolded region. Additionally, N- and C-terminal interactions stabilize HSP90 dimers, and molecular interactions in HSP90 dimers rearrange between the closed and semi-open states. Our findings provide valuable insight into the contributions of HSP90 and CDC37 in mediating client folding.

Comparative analysis of KRAS4a and KRAS4b splice variants reveals distinctive structural and functional properties.

KRAS, the most frequently mutated oncogene in human cancer, produces two isoforms, KRAS4a and KRAS4b, through alternative splicing. These isoforms differ in exon 4, which encodes the final 15 residues of the G-domain and hypervariable regions (HVRs), vital for trafficking and membrane localization. While KRAS4b has been extensively studied, KRAS4a has been largely overlooked. Our multidisciplinary study compared the structural and functional characteristics of KRAS4a and KRAS4b, revealing distinct structural properties and thermal stability. Position 151 influences KRAS4a's thermal stability, while position 153 affects binding to RAF1 CRD protein. Nuclear magnetic resonance analysis identified localized structural differences near sequence variations and provided a solution-state conformational ensemble. Notably, KRAS4a exhibits substantial transcript abundance in bile ducts, liver, and stomach, with transcript levels approaching KRAS4b in the colon and rectum. Functional disparities were observed in full-length KRAS variants, highlighting the impact of HVR variations on interaction with trafficking proteins and downstream effectors like RAF and PI3K within cells.

Allosteric Regulation of Switch-II Domain Controls KRAS Oncogenicity.

RAS proteins are GTPases that regulate a wide range of cellular processes. RAS activity is dependent on its nucleotide-binding status, which is modulated by guanine nucleotide exchange factors (GEF) and GTPase-activating proteins (GAP). KRAS can be acetylated at lysine 104 (K104), and an acetylation-mimetic mutation of K104 to glutamine (K104Q) attenuates the in vitro-transforming capacity of oncogenic KRAS by interrupting GEF-induced nucleotide exchange. To assess the effect of this mutation in vivo, we used CRISPR-Cas9 to generate mouse models carrying the K104Q point mutation in wild-type and conditional KrasLSL-G12D alleles. Homozygous animals for K104Q were viable, fertile, and arose at the expected Mendelian frequency, indicating that K104Q is not a complete loss-of-function mutation. Consistent with our previous findings from in vitro studies, however, the oncogenic activity of KRASG12D was significantly attenuated by mutation at K104. Biochemical and structural analysis indicated that the G12D and K104Q mutations cooperate to suppress GEF-mediated nucleotide exchange, explaining the preferential effect of K104Q on oncogenic KRAS. Furthermore, K104 functioned in an allosteric network with M72, R73, and G75 on the α2 helix of the switch-II region. Intriguingly, point mutation of glycine 75 to alanine (G75A) also showed a strong negative regulatory effect on KRASG12D. These data demonstrate that lysine at position 104 is critical for the full oncogenic activity of mutant KRAS and suggest that modulating the sites in its allosteric network may provide a unique therapeutic approach in cancers expressing mutant KRAS. SIGNIFICANCE: An allosteric network formed by interaction between lysine 104 and residues in the switch-II domain is required for KRAS oncogenicity, which could be exploited for developing inhibitors of the activated oncoprotein.

RAS-dependent RAF-MAPK hyperactivation by pathogenic RIT1 is a therapeutic target in Noonan syndrome-associated cardiac hypertrophy.

RIT1 is a RAS guanosine triphosphatase (GTPase) that regulates different aspects of signal transduction and is mutated in lung cancer, leukemia, and in the germline of individuals with Noonan syndrome. Pathogenic RIT1 proteins promote mitogen-activated protein kinase (MAPK) hyperactivation; however, this mechanism remains poorly understood. Here, we show that RAF kinases are direct effectors of membrane-bound mutant RIT1 necessary for MAPK activation. We identify critical residues in RIT1 that facilitate interaction with membrane lipids and show that these are necessary for association with RAF kinases and MAPK activation. Although mutant RIT1 binds to RAF kinases directly, it fails to activate MAPK signaling in the absence of classical RAS proteins. Consistent with aberrant RAF/MAPK activation as a driver of disease, we show that pathway inhibition alleviates cardiac hypertrophy in a mouse model of RIT1 mutant Noonan syndrome. These data shed light on the function of pathogenic RIT1 and identify avenues for therapeutic intervention.

Reduced dynamic complexity allows structure elucidation of an excited state of KRASG13D.

Localized dynamics of RAS, including regions distal to the nucleotide-binding site, is of high interest for elucidating the mechanisms by which RAS proteins interact with effectors and regulators and for designing inhibitors. Among several oncogenic mutants, methyl relaxation dispersion experiments reveal highly synchronized conformational dynamics in the active (GMPPNP-bound) KRASG13D, which suggests an exchange between two conformational states in solution. Methyl and 31P NMR spectra of active KRASG13D in solution confirm a two-state ensemble interconverting on the millisecond timescale, with a major Pγ atom peak corresponding to the dominant State 1 conformation and a secondary peak indicating an intermediate state different from the known State 2 conformation recognized by RAS effectors. High-resolution crystal structures of active KRASG13D and KRASG13D-RAF1 RBD complex provide snapshots of the State 1 and 2 conformations, respectively. We use residual dipolar couplings to solve and cross-validate the structure of the intermediate state of active KRASG13D, showing a conformation distinct from those of States 1 and 2 outside the known flexible switch regions. The dynamic coupling between the conformational exchange in the effector lobe and the breathing motion in the allosteric lobe is further validated by a secondary mutation in the allosteric lobe, which affects the conformational population equilibrium.

The ribosomal S6 kinase 2 (RSK2)-SPRED2 complex regulates the phosphorylation of RSK substrates and MAPK signaling.

Sprouty-related EVH-1 domain-containing (SPRED) proteins are a family of proteins that negatively regulate the RAS-Mitogen-Activated Protein Kinase (MAPK) pathway, which is involved in the regulation of the mitogenic response and cell proliferation. However, the mechanism by which these proteins affect RAS-MAPK signaling has not been elucidated. Patients with mutations in SPRED give rise to unique disease phenotypes; thus, we hypothesized that distinct interactions across SPRED proteins may account for alternative nodes of regulation. To characterize the SPRED interactome and evaluate how members of the SPRED family function through unique binding partners, we performed affinity purification mass spectrometry. We identified 90-kDa ribosomal S6 kinase 2 (RSK2) as a specific interactor of SPRED2 but not SPRED1 or SPRED3. We identified that the N-terminal kinase domain of RSK2 mediates the interaction between amino acids 123 to 201 of SPRED2. Using X-ray crystallography, we determined the structure of the SPRED2-RSK2 complex and identified the SPRED2 motif, F145A, as critical for interaction. We found that the formation of this interaction is regulated by MAPK signaling events. We also find that this interaction between SPRED2 and RSK2 has functional consequences, whereby the knockdown of SPRED2 resulted in increased phosphorylation of RSK substrates, YB1 and CREB. Furthermore, SPRED2 knockdown hindered phospho-RSK membrane and nuclear subcellular localization. We report that disruption of the SPRED2-RSK complex has effects on RAS-MAPK signaling dynamics. Our analysis reveals that members of the SPRED family have unique protein binding partners and describes the molecular and functional determinants of SPRED2-RSK2 complex dynamics.

Structural insights into the role of SHOC2-MRAS-PP1C complex in RAF activation.

RAF activation is a key step for signalling through the mitogen-activated protein kinase (MAPK) pathway. The SHOC2 protein, along with MRAS and PP1C, forms a high affinity, heterotrimeric holoenzyme that activates RAF kinases by dephosphorylating a specific phosphoserine. Recently, our research, along with that of three other teams, has uncovered valuable structural and functional insights into the SHOC2-MRAS-PP1C (SMP) holoenzyme complex. In this structural snapshot, we review SMP complex assembly, the dependency on the bound-nucleotide state of MRAS, the substitution of MRAS by the canonical RAS proteins and the roles of SHOC2 and MRAS on PP1C activity and specificity. Furthermore, we discuss the effect of several RASopathy mutations identified within the SMP complex and explore potential therapeutic approaches for targeting the SMP complex in RAS/RAF-driven cancers and RASopathies.

Consensus on the RAS dimerization hypothesis: Strong evidence for lipid-mediated clustering but not for G-domain-mediated interactions.

One of the open questions in RAS biology is the existence of RAS dimers and their role in RAF dimerization and activation. The idea of RAS dimers arose from the discovery that RAF kinases function as obligate dimers, which generated the hypothesis that RAF dimer formation might be nucleated by G-domain-mediated RAS dimerization. Here, we review the evidence for RAS dimerization and describe a recent discussion among RAS researchers that led to a consensus that the clustering of two or more RAS proteins is not due to the stable association of G-domains but, instead, is a consequence of RAS C-terminal membrane anchors and the membrane phospholipids with which they interact.

Structure of the SHOC2-MRAS-PP1C complex provides insights into RAF activation and Noonan syndrome.

SHOC2 acts as a strong synthetic lethal interactor with MEK inhibitors in multiple KRAS cancer cell lines. SHOC2 forms a heterotrimeric complex with MRAS and PP1C that is essential for regulating RAF and MAPK-pathway activation by dephosphorylating a specific phosphoserine on RAF kinases. Here we present the high-resolution crystal structure of the SHOC2-MRAS-PP1C (SMP) complex and apo-SHOC2. Our structures reveal that SHOC2, MRAS, and PP1C form a stable ternary complex in which all three proteins synergistically interact with each other. Our results show that dephosphorylation of RAF substrates by PP1C is enhanced upon interacting with SHOC2 and MRAS. The SMP complex forms only when MRAS is in an active state and is dependent on SHOC2 functioning as a scaffolding protein in the complex by bringing PP1C and MRAS together. Our results provide structural insights into the role of the SMP complex in RAF activation and how mutations found in Noonan syndrome enhance complex formation, and reveal new avenues for therapeutic interventions.

Exploring CRD mobility during RAS/RAF engagement at the membrane.

During the activation of mitogen-activated protein kinase (MAPK) signaling, the RAS-binding domain (RBD) and cysteine-rich domain (CRD) of RAF bind to active RAS at the plasma membrane. The orientation of RAS at the membrane may be critical for formation of the RAS-RBDCRD complex and subsequent signaling. To explore how RAS membrane orientation relates to the protein dynamics within the RAS-RBDCRD complex, we perform multiscale coarse-grained and all-atom molecular dynamics (MD) simulations of KRAS4b bound to the RBD and CRD domains of RAF-1, both in solution and anchored to a model plasma membrane. Solution MD simulations describe dynamic KRAS4b-CRD conformations, suggesting that the CRD has sufficient flexibility in this environment to substantially change its binding interface with KRAS4b. In contrast, when the ternary complex is anchored to the membrane, the mobility of the CRD relative to KRAS4b is restricted, resulting in fewer distinct KRAS4b-CRD conformations. These simulations implicate membrane orientations of the ternary complex that are consistent with NMR measurements. While a crystal structure-like conformation is observed in both solution and membrane simulations, a particular intermolecular rearrangement of the ternary complex is observed only when it is anchored to the membrane. This configuration emerges when the CRD hydrophobic loops are inserted into the membrane and helices α3-5 of KRAS4b are solvent exposed. This membrane-specific configuration is stabilized by KRAS4b-CRD contacts that are not observed in the crystal structure. These results suggest modulatory interplay between the CRD and plasma membrane that correlate with RAS/RAF complex structure and dynamics, and potentially influence subsequent steps in the activation of MAPK signaling.

Asynchronous Reciprocal Coupling of Martini 2.2 Coarse-Grained and CHARMM36 All-Atom Simulations in an Automated Multiscale Framework.

The appeal of multiscale modeling approaches is predicated on the promise of combinatorial synergy. However, this promise can only be realized when distinct scales are combined with reciprocal consistency. Here, we consider multiscale molecular dynamics (MD) simulations that combine the accuracy and macromolecular flexibility accessible to fixed-charge all-atom (AA) representations with the sampling speed accessible to reductive, coarse-grained (CG) representations. AA-to-CG conversions are relatively straightforward because deterministic routines with unique outcomes are achievable. Conversely, CG-to-AA conversions have many solutions due to a surge in the number of degrees of freedom. While automated tools for biomolecular CG-to-AA transformation exist, we find that one popular option, called Backward, is prone to stochastic failure and the AA models that it does generate frequently have compromised protein structure and incorrect stereochemistry. Although these shortcomings can likely be circumvented by human intervention in isolated instances, automated multiscale coupling requires reliable and robust scale conversion. Here, we detail an extension to Multiscale Machine-learned Modeling Infrastructure (MuMMI), including an improved CG-to-AA conversion tool called sinceCG. This tool is reliable (∼98% weakly correlated repeat success rate), automatable (no unrecoverable hangs), and yields AA models that generally preserve protein secondary structure and maintain correct stereochemistry. We describe how the MuMMI framework identifies CG system configurations of interest, converts them to AA representations, and simulates them at the AA scale while on-the-fly analyses provide feedback to update CG parameters. Application to systems containing the peripheral membrane protein RAS and proximal components of RAF kinase on complex eight-component lipid bilayers with ∼1.5 million atoms is discussed in the context of MuMMI.

Insights into the Cross Talk between Effector and Allosteric Lobes of KRAS from Methyl Conformational Dynamics.

KRAS is the most frequently mutated RAS protein in cancer patients, and it is estimated that about 20% of the cancer patients in the United States carried mutant RAS proteins. To accelerate therapeutic development, structures and dynamics of RAS proteins had been extensively studied by various biophysical techniques for decades. Although 31P NMR studies revealed population equilibrium of the two major states in the active GMPPNP-bound form, more complex conformational dynamics in RAS proteins and oncogenic mutants subtly modulate the interactions with their downstream effectors. We established a set of customized NMR relaxation dispersion techniques to efficiently and systematically examine the ms-µs conformational dynamics of RAS proteins. This method allowed us to observe varying synchronized motions that connect the effector and allosteric lobes in KRAS. We demonstrated the role of conformational dynamics of KRAS in controlling its interaction with the Ras-binding domain of the downstream effector RAF1, the first kinase in the MAPK pathway. This allows one to explain, as well as to predict, the altered binding affinities of various KRAS mutants, which was neither previously reported nor apparent from the structural perspective.

Machine learning-driven multiscale modeling reveals lipid-dependent dynamics of RAS signaling proteins.

RAS is a signaling protein associated with the cell membrane that is mutated in up to 30% of human cancers. RAS signaling has been proposed to be regulated by dynamic heterogeneity of the cell membrane. Investigating such a mechanism requires near-atomistic detail at macroscopic temporal and spatial scales, which is not possible with conventional computational or experimental techniques. We demonstrate here a multiscale simulation infrastructure that uses machine learning to create a scale-bridging ensemble of over 100,000 simulations of active wild-type KRAS on a complex, asymmetric membrane. Initialized and validated with experimental data (including a new structure of active wild-type KRAS), these simulations represent a substantial advance in the ability to characterize RAS-membrane biology. We report distinctive patterns of local lipid composition that correlate with interfacially promiscuous RAS multimerization. These lipid fingerprints are coupled to RAS dynamics, predicted to influence effector binding, and therefore may be a mechanism for regulating cell signaling cascades.

A Structure is Worth a Thousand Words: New Insights for RAS and RAF Regulation.

The RAS GTPases are frequently mutated in human cancer, with KRAS being the predominant tumor driver. For many years, it has been known that the structure and function of RAS are integrally linked, as structural changes induced by GTP binding or mutational events determine the ability of RAS to interact with regulators and effectors. Recently, a wealth of information has emerged from structures of specific KRAS mutants and from structures of multiprotein complexes containing RAS and/or RAF, an essential effector of RAS. These structures provide key insights regarding RAS and RAF regulation as well as promising new strategies for therapeutic intervention. SIGNIFICANCE: The RAS GTPases are major drivers of tumorigenesis, and for RAS proteins to exert their full oncogenic potential, they must interact with the RAF kinases to initiate ERK cascade signaling. Although binding to RAS is typically a prerequisite for RAF to become an activated kinase, determining the molecular mechanisms by which this interaction results in RAF activation has been a challenging task. A major advance in understanding this process and RAF regulation has come from recent structural studies of various RAS and RAF multiprotein signaling complexes, revealing new avenues for drug discovery.

The tumor suppressor activity of DLC1 requires the interaction of its START domain with Phosphatidylserine, PLCD1, and Caveolin-1.

BACKGROUND: DLC1, a tumor suppressor gene that is downregulated in many cancer types by genetic and nongenetic mechanisms, encodes a protein whose RhoGAP and scaffolding activities contribute to its tumor suppressor functions. The role of the DLC1 START (StAR-related lipid transfer; DLC1-START) domain, other than its binding to Caveolin-1, is poorly understood. In other START domains, a key function is that they bind lipids, but the putative lipid ligand for DLC1-START is unknown. METHODS: Lipid overlay assays and Phosphatidylserine (PS)-pull down assays confirmed the binding of DLC1-START to PS. Co-immunoprecipitation studies demonstrated the interaction between DLC1-START and Phospholipase C delta 1 (PLCD1) or Caveolin-1, and the contribution of PS to those interactions. Rho-GTP, cell proliferation, cell migration, and/or anchorage-independent growth assays were used to investigate the contribution of PS and PLCD1, or the implications of TCGA cancer-associated DLC1-START mutants, to DLC1 functions. Co-immunoprecipitations and PS-pull down assays were used to investigate the molecular mechanisms underlying the impaired functions of DLC1-START mutants. A structural model of DLC1-START was also built to better understand the structural implications of the cancer-associated mutations in DLC1-START. RESULTS: We identified PS as the lipid ligand for DLC1-START and determined that DLC1-START also binds PLCD1 protein in addition to Caveolin-1. PS binding contributes to the interaction of DLC1 with Caveolin-1 and with PLCD1. The importance of these activities for tumorigenesis is supported by our analysis of 7 cancer-associated DLC1-START mutants, each of which has reduced tumor suppressor function but retains wildtype RhoGAP activity. Our structural model of DLC1-START indicates the mutants perturb different elements within the structure, which is correlated with our experimental findings that the mutants are heterogenous with regard to the deficiency of their binding properties. Some have reduced PS binding, others reduced PLCD1 and Caveolin-1 binding, and others are deficient for all of these properties. CONCLUSION: These observations highlight the importance of DLC1-START for the tumor suppressor function of DLC1 that is RhoGAP-independent. They also expand the versatility of START domains, as DLC1-START is the first found to bind PS, which promotes the binding to other proteins.

RAS interaction with Sin1 is dispensable for mTORC2 assembly and activity.

RAS proteins are molecular switches that interact with effector proteins when bound to guanosine triphosphate, stimulating downstream signaling in response to multiple stimuli. Although several canonical downstream effectors have been extensively studied and tested as potential targets for RAS-driven cancers, many of these remain poorly characterized. In this study, we undertook a biochemical and structural approach to further study the role of Sin1 as a RAS effector. Sin1 interacted predominantly with KRAS isoform 4A in cells through an atypical RAS-binding domain that we have characterized by X-ray crystallography. Despite the essential role of Sin1 in the assembly and activity of mTORC2, we find that the interaction with RAS is not required for these functions. Cells and mice expressing a mutant of Sin1 that is unable to bind RAS are proficient for activation and assembly of mTORC2. Our results suggest that Sin1 is a bona fide RAS effector that regulates downstream signaling in an mTORC2-independent manner.

KRAS interaction with RAF1 RAS-binding domain and cysteine-rich domain provides insights into RAS-mediated RAF activation.

The first step of RAF activation involves binding to active RAS, resulting in the recruitment of RAF to the plasma membrane. To understand the molecular details of RAS-RAF interaction, we present crystal structures of wild-type and oncogenic mutants of KRAS complexed with the RAS-binding domain (RBD) and the membrane-interacting cysteine-rich domain (CRD) from the N-terminal regulatory region of RAF1. Our structures reveal that RBD and CRD interact with each other to form one structural entity in which both RBD and CRD interact extensively with KRAS. Mutations at the KRAS-CRD interface result in a significant reduction in RAF1 activation despite only a modest decrease in binding affinity. Combining our structures and published data, we provide a model of RAS-RAF complexation at the membrane, and molecular insights into RAS-RAF interaction during the process of RAS-mediated RAF activation.

Uncovering a membrane-distal conformation of KRAS available to recruit RAF to the plasma membrane.

The small GTPase KRAS is localized at the plasma membrane where it functions as a molecular switch, coupling extracellular growth factor stimulation to intracellular signaling networks. In this process, KRAS recruits effectors, such as RAF kinase, to the plasma membrane where they are activated by a series of complex molecular steps. Defining the membrane-bound state of KRAS is fundamental to understanding the activation of RAF kinase and in evaluating novel therapeutic opportunities for the inhibition of oncogenic KRAS-mediated signaling. We combined multiple biophysical measurements and computational methodologies to generate a consensus model for authentically processed, membrane-anchored KRAS. In contrast to the two membrane-proximal conformations previously reported, we identify a third significantly populated state using a combination of neutron reflectivity, fast photochemical oxidation of proteins (FPOP), and NMR. In this highly populated state, which we refer to as "membrane-distal" and estimate to comprise ∼90% of the ensemble, the G-domain does not directly contact the membrane but is tethered via its C-terminal hypervariable region and carboxymethylated farnesyl moiety, as shown by FPOP. Subsequent interaction of the RAF1 RAS binding domain with KRAS does not significantly change G-domain configurations on the membrane but affects their relative populations. Overall, our results are consistent with a directional fly-casting mechanism for KRAS, in which the membrane-distal state of the G-domain can effectively recruit RAF kinase from the cytoplasm for activation at the membrane.